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The interweaving chiral spirals (ICS), that is defined as superposition 
of differently oriented chiral spirals, is important for qualitative under- 
standings of the intermediate quark density region as well as quantitative 
estimates of the Quarkyonic region. We discuss how to construct the ICS, 
taking the (2+1) dimensional Fermi system as an example. We postu- 
late that the presence of the ICS would delay the occurrence of the chiral 
restoration as well as deconfinement phase transition, by tempering the 
growth of quark fluctuations. 

PACS numbers: 12:38Aw, 11.30.Rd 

1. Central Ideas: Qualitative Impacts of the ICS 

Recently, it has been argued that there is a new state of QCD matter 
at high baryon density and low to intermediate temperatures [1] (Figjl]). 

This novel state is called Quarkyonic matter, distinguished from nuclear 
matter by its bulk quantities such as pressure. The Fermi sea is mainly 
composed of quarks, not nucleons, in the region a little bit above fig ~ 
M]\f/Nc ~ Aqcd- This is because after the emergence of nucleons, a small 
change in rapidly enhances nucleon density, making strong short distance 
interactions among nucleons crucial. Thereby nucleons are not appropriate 
degrees of freedom to describe the bulk part of the Fermi sea - the quark 
picture is absolutely necessary to describe Quarkyonic matter. 

Quarkyonic matter should be also distinguished from conventional de- 
confined quark matter by its thermal and Fermi surface excitations. The 
excitations are confined, even after quarks are released from nucleons. A 
proper understanding of confined excitations is a basic starting point for any 
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Fig. 1. A speculated QCD phase diagram. 

discussions of phase structures, transport phenomena, and for construction 
of the effective Lagrangian. 

One of the most relevant observations in Ref. [1] is that the scale for 
the formation of the quark Fermi sea, /ig ~ Aqcd, and that for the decon- 

finement of the excitations, fiq ~ Nc^^'^ ^''Aqcd {d- spatial dimension), are 
conceptually different. The latter may be estimated by comparing quantum 
fluctuations of gluons with those of quarks near the Fermi surface. This 
observation was tested for (1+1) dimensional QCD, and it was argued that 
excitations are always confineqj indendently of Nc, while the pressure is 
saturated by free quark contributions [2]. 

Now let us ask what happens to chiral symmetry. This issue is very 
important for qualitative understandings of the intermediate density region 
as well as quantitative estimates for the Quarkyonic region. Below we will 
argue that chiral symmetry is spontaneously broken by inhomogeneous chiral 
spirals, tempering the growth of the quark fluctuations at finite density. 

Sometimes it is said that Quarkyonic matter is defined as a chiral sym- 
metric confined matter, although chiral symmetry was not a primary is- 
sue in the original proposal of Ref. |T]. If one sticks to this definition, it 
might be very difficult to imagine the existence of Quarkyonic matter at 
A'^c = 3. Indeed, once the quarks near the Fermi surface become gapless, 
it would largely enhance quark fiuctuations in addition to the phase space 
enhancement at finite density. Recent results of the functional renormal- 
ization group application to the PNJL model for not very high density [3] 
presumablly should be interpreted in this context. 

A general tendency of model analyses, without depending on whether 
models are confining or not, suggests that the chiral restoration occurs 



^ In spatial one dimension, the phase space for quark fluctuations is always the same 
as vacuum case, so the confining force is never modified. 
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Fig. 2. Three types of chiral pairings. 

shortly after the formation of the quark Fermi sea, or ~ Aqcd [II- This 
trend may be understood by observing that the creation of anti-quarks, 
which are ingredients of the usual chiral condensate, costs more energy for 
the larger quark Fermi sea. This is so because the particle in the Dirac sea 
must go above the Fermi surface to avoid the Pauli-blocking (Figj2]). 

At finite density, however, more proper ingredients of condensates are 
particle-holes near the Fermi surface. The excitations near the Fermi surface 
naturally have large momenta, \f)\ ~ /Xg, but they do not cost additional 
kinetic energy much, compared to the energy before excitations. Shown in 
Figj2]are the exciton and density wave pairings. 

In an exciton case, the total momenum of a pair is ~ 0, but the relative 
momentum between a particle and a hole is ~ 2^g. In a confining model, 
such a pairing accompanies a large string, costing large potential energy. 

In a density wave case, while the total momenum of a pair is ~ 2/ig, 
a particle and a hole co-move without forming a large string. Thus the 
density wave pairing should be energetically favored compared to the exciton 
pairing. (Arguments based on confining model picture is useful but not 
indispensable, though. See below.) 

Actually the chiral density wave solution can be always interpreted as 
the chiral spirals. A key observation is that once we have a condensation 
of a pair moving to, say, +z-direction, there is also a pair moving to —z 
direction. Mathmatically, one can project out fermion components moving 
to ±2; directions by operating the projection matrices [5], 

_ 1 ±707^ 

V'i = — ^ — ^ ■ (1) 
Then we have two types of the chiral condensates, 

~ Ae^i^'^^ , (V^+V'-) ~ Ae-2iM,^ ^ (2) 
whose sum and difference give 



{tpij)) ~ A cos(2/ig2;) , (V'iToTzV') ~ A sin(2^gz) , (3) 
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with a fixed radius of A of the order Aq(-.j-). These condensates obviously 
break the chiral symmetry, translational invariance, rotational invariance, 
and the second condensate further breaks parity locally. 

Here once again we emphasize that in the presence of the Fermi sea, 
large momenta naturally appear without costing much excitation energies, 
so one should not be surprised at the emergence of condensations with 
relatively large momenta. Indeed, analyses of both non-confining [3 [71 [HI [9] 
and confining models [5j suggest that the chiral spiral solution overtakes the 
homogeneous solution. 

Rather a more nontrivial question is related to the fact that the chiral 
spiral must have a particular orientation. Let us ask: Can chiral pairs be 
formed in such a way to cover the entire Fermi surface, and can differently 
oriented chiral spirals be interweaved in a consistent way? 

To answer to these questions will be very important for considerations 
about whether the chiral symmetry breaking may survive after taking the 
color superconductivity into account. If only single chiral spiral in one 
particular direction were possible, we could employ a less number of pairs 
for chiral condensations than we do for diquark condensations. If this would 
be the case, the color superconducting phase would overtake a single chiral 
spiral [?]• But instead we shall suggest a possibility that the ICS appears 
as far as the nonperturbative gluon exchange survives. 

A possibility of the ICS has been qualitatively discussed in a confining 
model LlOj, although whole aspects about the ICS were not fully explored 
because of the technical difficulties in treating the deep infrared region of 
the gluon exchange. Actually, however, key aspects about the ICS may be 
extracted without using the confinement, and, in fact, are rather robust to 
the detailed behaviors in the deep infrared region. 

Below we shortly highlight how to construct the ICS for the non-confining 
model in (2-1-1) dimensions, together with the parametric estimates of sev- 
eral effects. How to handle the corrections, relations with the previous works 
[Sl[n], etc., have been comprehensively discussed in a recent paper p!T], so 
an interested reader should consult it for details. 



2. How to Construct the ICS in (2-(-l) dimensions 

As an example of the ICS, we consider a (2-1-1) dimensional Fermi sccj^ 
at T = 0, purely because possible shapes of the Fermi surface are relatively 
simple. An extension of our treatments to higher dimensional systems are 
technically nontrivial but conceptually straightforward. 

^ Precisely speaking, in odd dimensional space-time, the chirality is not defined, so 
the terminology "chiral" spirals may be a little bit misleading. But the mechnism to 
generate spirals do not depend on this fact. 
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Particle-hole combinations 
for one patch chiral spirals 



Fig. 3. A 3-patches = 6-wedges case. 

We first divide tlie Fermi sea into 2N^ wedges (FigjsJ). A wegde in one 
side of the Fermi sea and a wedge in the opposite side are regarded as a 
pair, and we call it one patch domain. We denote the height of each wedge 
as Q, and its open angle as 20 = n/Np. Each patch will generate a single 
chiral spiral. The variables Q and Q are variational parameters which will 
be optimized in such a way to minimize the total free energy. 

Below we will use the canonical emsemble with specifying quark number 
density or pp, because the presentation of ideas becomes simpler than the 
grand canonical case. Then the Fermi volume conservation can be used to 
rewrite Q as a function ol pp and 0. In the canonical emsemble, our goal 
is to minimize the total energjj^by choosing the optimal value of 0. 

Essentially, the total energy and the shape of the Fermi surface will be 
determined by balancing the following energy costs and gains (FigQ: 

• (i) The kinetic energy cost arising from the deformation of the Fermi 
surface, which weakly depends upon the condensation effects. This 
contibution becomes dominant for large 0. 

• (ii) The energy gain in a single patch from the condensation effects. 
The condensation effects bend down the single particle dispersion. 
Then particles occupy smaller energy orbits, reducing the total single 
particle energy. In the following the quark mass gap will be denoted 
as M. This contribution is less sensitive to the angle compared to 
other contributions. 

• (iii) Coherent interactions among differently oriented chiral spirals (we 
call it inter-patch interaction) which cost energy. They originate from 
the condensation effects, so become less important as condensates get 
smaller. It becomes increasingly important for small due to the 
enhanced number of inter-patch interactions. 



At T = 0, the minimization of tiie free energy in the grand canonical emsemble is 
thermodynamically equivalent to the minimization of the total energy in the canonical 
emsemble. 
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■ Cost : Deformation 

(dominant for large 0) 




equal volume 



■ Gain : Mass gap origin 




Condensation effects 




■ Cost : Interactions among CSs 



(dominant for small 0) 



Condensate - Condensate int. 



Fig. 4. The energy gains and costs. 



Below we shall give estimates for these contributions, starting with the 
thereotically clean set up, that is, the leading order (LO) of the ex- 
pansion and the high density expansion in powers of Aqcd/pf- Although 
LO results are not directly applied to the phenomenologically interesting 
region, it is not difficult to specify which effects will grow at lower density, 
and how results will be modified qualitatively. Indeed, many features of the 
previous works [8l [9] , which have been numerically done for relatively lower 
density, can be understood from this analytic framework. 

For explicit estimates, it is indispensable to introduce models. In par- 
ticular, the momentum dependence of the interaction strongly affects the 
estimate of the inter-patch interactions. We will explain the features of the 
model, then move to parametric estimates of several effects. 



In many studies of the intermediate density, typically the 4-Fermi in- 
teraction has been used, with the ultraviolet (UV) cutoff on the quark 
momenta. But this approach will be problematic when the quark Fermi 
momentum pp becomes closer to the UV cutoff. Thus we have to use al- 
ternative descriptions at the intermediate density region. Our model is the 
following non-local 4-Fermi interaction, 



2.1. A model 




(4) 



where Op^k = ^(^qcd ~ iP~ ^)^)) momentum integration is for q,p,k. 
The large- A'^c QCD is mimicked as follows. The one-gluon exchange includ- 
ing non-perturbative effects are shown in Fig. ^a). Its strength damps as 
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the momentum transfer becomes large. We roughly take into account this 
property by introducing a step function, 6{Aqqy) — (p — k)"^)^ keeping the 
interaction strength constant. 

In Fig. [sjl^b), we show the color line representation to illustrate how the 
one-gluon exchange interaction should be contracted into a four- Fermi type 
interaction. Taking into account features in Figs. [Sj^a) and[5]^b), we arrive 
at a simple model described in Eq. Q and Fig. [sj^c). 

The main consequence of our form factor treatments can be best seen 
in the Schwinger-Dyson or gap equations (See Fig. [5]). We illustrate it for 
zero density case. After picking up a residue, we arrive at 

Mm = /(l^ ^ ^P.^ • = ^M(kY + P^ (5) 

Note that because of Op ^^ the contributions to the mass gap M{p) comes 
from the integral around p. Putting in a different way, a particle with p is af- 
fected by the condensate made of particles and anti-particles with momenta 
close to p. So quark-condensate interactions are local in momentum space. 
Also it is clear that at very large \p\, the integrand quickly damps, leading 
to small M{pj. The chiral restoration occurs for high energy excitations. 

By applying this argument, one can derive several conclusions about 
finite density. The low energy excitations appear near the Fermi surface, 
so that the chiral symmetry is violated near the Fermi surface, but it is 
gradually restored in the region far from the Fermi surface. 

Thanks to the locality in momentum space, the inter-patch interactions 
among condensates occur only if their momentum domains are close. In the 
construction of the ICS, this property is essential to restrict interactions 
among differently oriented chiral spirals only near the patch boundaries. 
If this locality is absent, inter-patch interactions would occur everywhere 
in the Fermi surface destroying the chiral spirals one another, and would 
reduce the quark mass gap considerably [9]. 

Now we are ready to give qualitative estimates of several contributions. 

2.2. The energy cost: Deformation energy 

The contribution (i) is rather easily estimated. In one wedge, Q is fixed 
by the Fermi volume conservation, 

ple = Q'^tane, (6) 

and the difference between the Deformed Fermi sea and the spherical Fermi 
sea can be calculated as 

AWm.(e) ~ iVc • p|e^ + 0(66) ^ 



8 



wroclaw'kojo printed on November 21, 2011 




Fig. 5. (a) The non-perturbative gluon exchange which is supposed to damp quickly 
in the UV region, (b) The color line representation of the one-gluon exchange, (c) 
Our effective four-Fermi interaction including form factor effects. 




small momentum transfer form factor 



Fig. 6. The leading self-energy diagram at zero density. (Left) The diagram in 
terms of QCD dynamics with the Rainbow Ladder approximation. (Right) The 
corresponding diagram in our model. 

where we have assumed <C 1 and M/pp <C 1. Interestingly, term 
disappears after the subtraction of the spherical Fermi sea. 



2.3. The energy gain: Single particle energy 

The solution of the gap equation gives the mass gap, M ~ Aqcd, opened 
near the Fermi surface within distance of ~ Aqcd- On the other hand, 
quarks outside of this domain are not strongly affected by the condensates. 
Thus the energy reduction after summing up all patch contributions is 

Afcond.(0) ~ A^c • (-M) X (Aqcd -PF tan 9) x iVp ~ -iY^A^cDPF , (8) 

where —M and Aqcd ■ p_Ftaii20 characterizes the bending down of single 
particle energy, a number of particles acquiring the mass gap within a single 
patch, respectively. A sum of all patch contributions is approximately 
independent, and does not play a relevant role in the optimization of 0. 
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2.4- The energy cost: Inter-patch interactions among condensates 

Since the single particle and the condensate interact only if their do- 
mains in momentum space are close one another, so interactions between 
differently oriented chiral spirals occur only near the patch boundaries. The 
strength of the interaction is proportional to M^. Taking into account 
the phase space where interactions occur, and counting a number of patch 
boudaries, we estimate the energetic cost as 

Afint. (6) - iVc • M^Aqcd X TVp ~ • (9) 

Note that the contributions are proportional to 1/G, so the creation of 
condensates are not favored for very small 0. Actually, near the patch 
boundaries the quark mass gap becomes effectively smaller. 



2.5. The optimimal value of Q 

Now we can optimize 6 by differentiating the total energy. Since the 
single particle contributions do not strongly depend on details of 0, the 
optimal value of G is essentially determined by balancing the deformation 
energy and the inter-patch interaction energy: 

^ (AWm. + ^£int.) ~ 4p|e3 - ^ ~ , (10) 

which determines the optimal value of 6 as 

Therefore, a number of chiral spirals, A?p, increases as density does. Since 
A'p is an integer, the phase transition occurs discontinuously. 

With this value of G, the deformation and interaction energies are ^ 
Pp^Aq(^^, smaller than the single particle contribution, ~ — pi^Aq^^,. 



3. Summary 

We have argued why the ICS is potentially relevant, and shown how to 
construct it by taking the (2-1-1) dimensional Fermi system as an example. 

The chiral spirals is a mechanism that can generate the quark mass gap 
of O(Aqcd) even after the formation of the quark Fermi sea. We expect that 
the mass gap would temper the growth of the quark fluctuations near the 
Fermi suface, shifting chiral restoration and deconfinement lines to higher 
temperature and density than those predicted under the assumption of the 
homogeneous condensates. 
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We have argued the ICS only at T = 0, but for phenomenological apph- 
cations to the RHIC low energy scan or future FAIR and NICA experiments, 
the extension of the present results to T 7^ are absolutely necessary. We 
expect that some qualitative changes occur at some temperature. Such an 
extension will be presented in near future. 
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